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ABSTRACT 

Hydrodynamic simulations of the merger of stellar mass black hole - neutron star binaries (BH/NS) are com- 
pared with mergers of binary neutron stars (NS/NS). The simulations are Newtonian, but take into account the 
emission and backreaction of gravitational waves. The use of a physical nuclear equation of state allows us to 
include the effects of neutrino emission. For low neutron star to black hole mass ratios the neutron star transfers 
mass to the black hole during a few cycles of orbital decay and subsequent widening before finally being disrupted, 
whereas for ratios near unity the neutron star is already distroyed during its first approach. A gas mass between 
~ 0.3Mq and ~ 0.7 Mq is left in an accretion torus around the black hole and radiates neutrinos at a luminosity of 
several lO^'' erg/s during an estimated accretion time scale of about 0.1 s. The emitted neutrinos and antineutrinos 
annihilate into pairs with efficiencies of 1-3% percent and rates of up to ~ 2 x lO^^erg/s, thus depositing an 
energy E^c^ ^10^' erg above the poles of the black hole in a region which contains less than 10"^ Mq of baryonic 
matter. This could allow for relativistic expansion with Lorentz factors around 100 and is sufficient to explain 
apparent burst luminosities ~ E^ij/ifat^) up to several lO^-'ergs"' for burst durations « 0.1-ls, if the 7 
emission is collimated in two moderately focussed jets in a fraction fa = 25VI/{Ati) k, 1/100 — 1/10 of the sky. 

Subject headings: binaries: close — black hole physics — gamma-rays: bursts — stars: neutron 



{vv) annihilation as potential source of energy for GRBs. 

2. NUMERICAL METHODS 

The three-dimensional hydrodynamic simulations were per- 
formed with a Eulerian PPM code using four levels of nested 
cartesian grids which ensure a good resolution near the center of 
leasi lor me suDciass or less complex anu less energetic snori m^^s and a large computational volume simultaneously. Each 
! or,H hn.r-A K,„-ctc t\/i^,^ Mo,-o„o„ Sr Dit-o^ \QQA\ ..Atu H.^-oti^^c grid had 64^ zoncs, the size of the smallest zone was 0.64 or 

0.78 km in case of NS/NS and 1 .25 or 1 .5 km for BH/NS merg- 
ers. The zone sizes of the next coarser grid levels were doubled 
to cover a volume of 328 or 400km side length for NS/NS and 
640 or 768km for BH/NS simulations. GW emission and its 
backreaction on the hydrodynamics were taken into account by 
the method of Blanc het, Damour, & Schafer (1990) (see Ruf- 
fert, Janka, & Schafer 1996). The neutrino emission and corre- 
sponding energy and lepton number changes of the matter were 
calculated with an elaborate neutrino leakage scheme (Ruffert, 
Janka, & Schafer 1996), and vv annihilation around the merger 
was evaluated in a post-processing step (Ruffert et al. 1997). 

3. SIMULATIONS 

Table 1 gives a list of computed NS/NS and BH/NS merger 
models. Besides the baryonic mass of the neutron star and the 
mass of the black hole, the spins of the neutron stars were var- 
ied. "Solid" means synchronously rotating stars, "none" irro- 
tational cases and "anti" counter-rotation, i.e., spin vectors op- 
posite to the direction of the orbital angular momentum. The 
cool neutron stars have a radius of about 15 km (Ruffert, Janka, 
& Schafer 1996) and the runs were started with a center-to- 
center distance of 42^6 km for NS/NS and with 47 km in case 
of BH/NS for Mbh = 2 .5 Mq , 57 km for Mbh = 5 Mq and 72 km 
forMfiH = IOMq. The simulations were stopped at a time fsim 
between 10ms and 20ms. The black hole was treated as a point 
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1. INTRODUCTION 

BH/NS and NS/NS mergers are discussed as promising 
candidates for the origin of gamma-ray bursts (GRBs) (e.g., 
BUnnikov et al. 1984; Eichler et al. 1989; Paczyhski 
1991; Narayan, Piran, & Shemi 1992; Meszaros 1999; Fryer, 
Woosley, & Hartmann 1999; Bethe & Brown 1998, 1999), at 
least for the subclass of less complex and less energetic short 
and hard bursts (Mao, Narayan, & Piran 1994) with durations 
of fractions of a second (Popham, Woosley, & Fryer 1999; Ruf- 
fert & Janka 1999). Optical counterparts and afterglows of this 
subclass have not yet been observed. Due to the presence of a 
region of very low baryon density above the poles of the black 
hole, BH/NS mergers are considered as more favorable sources 
than NS/NS mergers (e.g., Portegies Zwart 1998; Brown et al. 
1999). 

Previous Newtonian SPH simulations of BH/NS mergers us- 
ing a polytropic equation of state indicate that the neutron star 
may slowly lose gas in many mass transfer cycles (Kluzniak 
& Lee 1998; Lee & Kluzniak 1998,1999). Whether dynamical 
instability sets in at a minimum separation (Rasio & Shapiro 
1994; Lai, Rasio, & Shapiro 1994) or whether stable Roche 
lobe overflow takes place, however, can depend on the neutron 
star to black hole mass ratio (Bildsten & Cutler 1992) and the 
properties of the nuclear equation of state, expressed by the adi- 
abatic index (Uryu & Eriguchi 1999). 

In this Letter we report about the first Newtonian BH/NS 
merger simulations (Eberl 1998) which were done with a real- 
istic nuclear equation of state (Lattimer & Swesty 1991) and 
which therefore yield information about the thermodynamic 
evolution and the neutrino emission. They allow one to com- 
pare the strength of the gravitational wave (GW) emission rela- 
tive to NS/NS mergers and to investigate neutrino-antineutrino 
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mass at the center of a sphere with radius Rs = 2GMbh/c^ which 
gas could enter unhindered. Its mass and momentum were up- 
dated along with the accretion of matter. Model TNIO, which 
is added for comparison, is a continuation of the NS/NS merger 
model B64 where at time fsim = 10 ms the formation of a black 
hole was assumed and the accretion was followed for another 
5 ms until a steady state was reached (Ruffert & Janka 1999). 

4. RESULTS 

4.1. Evolution of BH/NS mergers 

Due to the emission of GWs the orbital separation decreases. 
During its first approach, the neutron star transfers matter to 
the black hole at huge rates of several 100 up to ^ IOOOMq/s. 
Within 2-3 ms it loses 50-75% of its initial mass. In case of the 
2.5 M0 black hole the evolution is catastrophic and the neutron 
star is immediately disrupted (Lattimer & Schramm 1974). A 
mass of 0.2-0.3 Mq remains in a thick disk around the black 
hole (Md in Table 2). In contrast, the orbital distance increases 
again for Mbh = 5Mq and Mbh = IOMq and a significantly 
less massive neutron star begins a second approach. Again, the 
black hole swallows gas at rates of more than IOOM0/S. Even 
a third cycle is possible (Fig. 1). Finally, at a distance dns and 
time fns the neutron star with a mass of M™" is destroyed and 
most of its mass ends up in an accretion disk (Table 2). (In 
case of NS/NS mergers ?ns means the time when the two den- 
sity maxima of the stars are one stellar radius, i.e., dns = 15 km, 
apart). 

The increase of the orbital separation is connected with a 
strong rise of the specific (orbital) angular momentum of the gas 
(Fig. 1). Partly this is due to the fact that the black hole can cap- 
ture gas with low specific angular momentum first, but mainly 
because only a fraction of the orbital angular momentum of the 
accreted gas is fed into spinning up the black hole. This frac- 
tion, which is lost for the orbital motion, is proportional to the 
quantity a in Fig. 2. Figure 2 is based on the parameterized 
analysis of non-conservative mass-transfer by Podsiadlowski, 
Joss, & Hsu (1992) (see also Fryer et al. 1999) assuming that 
mass ejection from the system is negligible. It shows that dis- 
regarding GW emission, the orbital separation can increase for 
small initial black hole mass only after the neutron star has lost 
much mass, while for larger initial Mbh and smaller a orbital 
widening is easier Without GWs the separation increases when 
a < (Mbh-Mns)/(Mbh+Mns). Including angular momentum 
loss by GWs in the point-mass approximation and using the 
mass-loss rates from the hydrodynamic models (dashed lines in 
Fig. 2) yields a quahtative understanding of the behavior visible 
in Fig. 1 and suggests that a is between 0.2 and 0.5. 

During the merging a gas mass AMej of ~ 10~'^ Mq (in case 
of counter-rotation and Mbh = 2.5 Mq) to ~ 0.1 M0 (corota- 
tion and Mbh = IOMq) is dynamically ejected (Table 2). In the 
latter case the associated angular momentum loss is about 7%, 
in all other cases it is less than 5% of the total initial angular 
momentum of the system. Another fraction of up to 24% of 
the initial angular momentum is carried away by GWs. In Ta- 
ble 2 the rotation parameter a = Jc/(GM^) is given for the initial 
state of the binary system (ai) and at the end of the simulation 
(af) for the renmant of NS/NS mergers or for the black hole in 
BH/NS systems, respectively, provided the black hole did not 
have any initial spin. When the whole disk mass Mj has been 
swallowed by the Kerr black hole, a final value (Table 2) 
will be reached in case of the accretion of a corotating, thin disk 
with maximum radiation efficiency. 



The phase of largest mass flow rate to the black hole (be- 
tween 2 and 5 ms after the start of the simulations) is connected 
with a maximum of the GW luminosity Lqw which reaches up 
to 7 X 10^^ erg/s (Table 1). The peak values of Lqw and the wave 
amplitude rh (for distance r from the source) increase with the 
black hole mass. The total energy £gw radiated in GWs can be 
as much as 0. 1 Mqc^ for Mbh = IOMq . 

4.2. Neutrino Emission and GRBs 

Compressional heating, shear due to numerical viscosity, and 
dissipation in shocks heat the gas during accretion to maximum 
temperatures kT"^'^ of several lOMeV. Average temperatures 
are between 5 and 20MeV, the higher values for the less mas- 
sive and more compact black holes. At these temperatures and 
at densities of 10'"-10'^g/cm^ in the accretion flow, electrons 
are non-degenerate and positrons abundant. Electron neutrinos 
and antineutrinos are therefore copiously created via reactions 
p+e~ ^ n+Ve andn+e^ ^ p+i'e sad dominate the neutrino en- 
ergy loss from the accreted matter Dense and hot neutron mat- 
ter is not completely transparent to neutrinos. By taking into 
account the finite diffusion time, the neutrino trapping scheme 
limits the loss of energy and lepton number. 

In Table 1 maximum and average values of the luminosi- 
ties (L™" and Lff., respectively, the latter in brackets) in the 
simulated time intervals are hsted for i/g and i?e and for the 
sum of all heavy-lepton neutrinos. The latter are denoted by 
Vx = v^j_,v^,Vr, Vr and are mainly produced by e^e~ annihi- 
lation. The total neutrino luminosities L^(f) (Fig. 3) fluctuate 
strongly with the varying mass transfer rate to the black hole 
during the cycles of orbital decay and widening (compare with 
Fig. 1). The total energy E^, radiated in neutrinos in 10-20ms is 
typically several 10^' erg. Time averages of the mean energies 
(e) of the emitted neutrinos are ~ 15MeV for 20MeV for v^, 
and 30MeV for v^. Luminosities as well as mean energies, in 
particular for smaller black holes, are significantly higher than 
in case of NS/NS mergers. 

At the end of the simulations, several of the BH/NS models 
have reached a steady state, characterized by only a slow growth 
of the black hole mass with a nearly constant accretion rate. 
Corresponding rates are given in Table 2 and are several 
Mq/s. From these we estimate torus Ufe times face = Mj/Md of 
50-1 50ms. Values with > and < signs indicate cases where the 
evolution and emission are still strongly time-dependent at fsim- 
In these cases the accretion torus around the black hole has also 
not yet developed axial symmetry. In all other cases the effec- 
tive disk viscosity parameter aeff ~ Vr/vKepier ~ 3-\/6/?s/(?accc), 
evaluated at a representative disk radius of 3/?s = 6GMbh/c^, 
has the same value, 4-5 x 10"'. This value is associated with 
the numerical viscosity of the hydro code (which solves the Eu- 
ler equations) for the chosen resolution. The further disk evo- 
lution is driven by the angular momentum transport mediated 
by viscous shear forces, which determines the accretion rate. 
The physical value of the disk viscosity is unknown. The nu- 
merical viscosity of our code, however, is in the range where 
the viscous energy dissipation and the energy emission by neu- 
trinos should be roughly equal, i.e., where the conversion ef- 
ficiency qu = (L,y)/(MdC^) of rest-mass energy to neutrinos is 
nearly maximal (see Ruffert et al. 1997; Ruffert & Janka 1999). 

Assuming that the average neutrino luminosity (L,^) at fsim is 
representative for the subsequent accretion phase, we obtain for 
qv numbers between 4 and 6% and total energies ~ (^i/)^acc 
around 3 x 10^^ erg (Table 2). Annihilation of neutrino pairs, 
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uV e^e~, deposits energy at rates up to fi^p ^ 2 x lO^^erg/s 
in the vicinity of the black hole (Fig. 4). This corresponds to 
total energies ^ Euvhcc as high as ~ 10^' erg and annihi- 
lation efficiencies = -E^p/ {L^) of 1-3%. These estimates 
should not change much if the different effects of general rela- 
tivity on vv annihilation are taken into account in combination 
(Ruffert & Janka 1999; Asano & Fukuyama 1999), but gen- 
eral relativistic simulations of the merging are very important. 
More energy could be pumped into the e^7 fireball when the 
black hole rotates rapidly (Popham, Woosley, & Fryer 1999) or 
if magnetic fields are able to tap the rotational energy of the ac- 



cretion torus and of the black hole with higher efficiency than 
vv annihilation does (Blandford & Znajek 1977). This seems to 
be necessary for the long and very energetic GRBs (Meszaros, 
Rees, & Wijers 1999; Brown et al. 1999; Lee, Wijers, & Brown 
1999). 
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TABLE 1 

GRAVITY WAVES AND NEUTRINOS FROM NS/NS AND BH/NS MERGING 



Model 


Type 


Masses 


Spin 




rmax 






J max(av) 


J max(av) 


J max{av) 

Si/, 
















Wo) 




(ms) (104^f)'' 


(iU cm 


) (foe) (100 if) (100 if) (100 ^f) (foe) (MeV) (MeV) (MeV) (MeV) 


S64 


NS/NS 




solid 


10 


0.7 


5.5 


14 


0.3(0.2) 


0.9(0.5) 


0.3(0.2) 


0.8 


35 


12 


18 


26 


D64 


NS/NS 


1.2+1.8 


solid 


13 


0.4 


5.5 


13 


0.5(0.3) 


1.3(0.8) 


0.7(0.4) 


1.1 


35 


13 


19 


27 


V64 


NS/NS 


1.6+1.6 


anti 


10 


1.2 


6.0 


23 


1.1(0.5) 


2.6(1.3) 


0.7(0.3) 


1.9 


69 


13 


19 


27 


A64 


NS/NS 


1.6+1.6 


none 


10 


2.1 


8.6 


52 


0.9(0.5) 


2.6(1.3) 


1.4(0.6) 


2.3 


39 


12 


18 


26 


B64 


NS/NS 


1.6+1.6 


solid 


10 


2.1 


8.9 


37 


0.6(0.4) 


1.8(1.1) 


0.9(0.4) 


1.8 


39 


13 


19 


27 


TNIO 


BH/AD 


2.9+0.26 


solid 


15 








0.5(0.4) 


1.3(0.9) 


0.6(0.2) 


0.8 


15 


9 


13 


21 


C2.5 


BH/NS 


2.5+1.6 


anti 


10 


2.3 


9.9 


32 


1.5(0.5) 


7.3(2.5) 


5.2(1.9) 


4.5 


74 


16 


22 


31 


A2.5 


BH/NS 


2.5+1.6 


none 


10 


2.0 


9.9 


50 


1.8(0.5) 


6.4(2.2) 


3.1(1.3) 


3.6 


65 


15 


22 


31 


B2.5 


BH/NS 


2.5+1.6 


solid 


10 


2.1 


9.6 


61 


0.9(0.3) 


6.5(1.7) 


3.6(0.9) 


2.5 


61 


14 


21 


29 


C5 


BH/NS 


5.0+1.6 


anti 


15 


3.9 


13.0 


50 


0.7(0.4) 


3.8(1.6) 


2.5(1.1) 


4.5 


46 


15 


20 


29 


A5 


BH/NS 


5.0+1.6 


none 


20 


3.2 


14.8 


102 


0.7(0.2) 


4.4(1.5) 


2.8(0.8) 


4.5 


51 


16 


24 


31 


B5 


BH/NS 


5.0+1.6 


solid 


15 


3.4 


14.5 


95 


0.6(0.2) 


3.7(1.1) 


2.5(0.6) 


2.9 


44 


14 


21 


28 


CIO 


BH/NS 


10.0+1.6 


anti 


10 


7.1 


21.9 


123 


0.4(0.1) 


2.5(0.4) 


1.2(0.1) 


0.6 


51 


14 


19 


24 


AlO 


BH/NS 


10.0+1.6 


none 


10 


6.9 


26.2 


168 


0.2(0.1) 


2.5(0.5) 


1.2(0.2) 


0.7 


50 


14 


20 


26 


BIO 


BH/NS 


10.0+1.6 


solid 


10 


7.3 


26.2 


163 


0.4(0.1) 


2.5(0.8) 


1.4(0.2) 


1.1 


52 


13 


18 


24 



1 foe =10 erg (fifty one erg). 



TABLE 2 

DISK FORMATION AND NEUTRINO ANNfflmATION 



Model 


^ns 






AMej 




Md 


^acc 


Q!eff 




at 


„oo 
"BH 












F - 




(ms) (km) (Mq) (Mq/100) (Mq) (Mq/s) 


(ms) 


(10-3) 








(100 ^f)" 


(foe/s) 


(%) 


(%) 


(foe) 


(foe) 


S64 


2.8 


15 




2.0 










0.98 


0.75 




1.5 


1 




1 






D64 


7.3 


15 




3.8 










0.87 


0.69 




2 


2 




1 






V64 


3.7 


15 




0.0085 










0.64 


0.49 




4 


9 




2 






A64 


1.7 


15 




0.23 










0.76 


0.55 




5 


9 




2 






B64 


1.6 


15 




2.4 










0.88 


0.63 




3 


7 




2 






TNIO 










0.26 


5 


53 


4 




0.42 


0.59 


1.2 


0.5 


1.3 


0.4 


7 


0.03 


C2.5 


2.6 


11 


0.78 


0.01 


0.26 


6 


43 


4 


0.65 


0.47 


0.60 


7 


20 


6 


3 


30 


0.9 


A2.5 


4.3 


18 


0.78 


0.03 


0.33 


< 14 


>24 


<8 


0.67 


0.39 


0.56 


7 


20 


>3 


3 


> 17 


>0.5 


B2.5 


6.0 


23 


0.78 


0.2 


0.45 


<35 


> 13 


< 14 


0.69 


0.38 


0.61 


7 


20 


> 1 


3 


>9 


>0.3 


C5 


9.1 


76 


0.40 


2.5 


0.38 


5 


76 


5 


0.44 


0.27 


0.42 


4 


8 


4 


2 


30 


0.6 


A5 


16.3 


65 


0.52 


2.5 


0.49 


6 


82 


4 


0.45 


0.17 


0.37 


4 


8 


4 


2 


33 


0.7 


B5 


10.8 


79 


0.50 


5.6 


0.45 


6 


75 


5 


0.46 


0.19 


0.38 


4 


8 


4 


2 


30 


0.6 


CIO 


8.0 


96 


0.65 


2.2 


0.67 


< 10 


>67 


< 11 


0.24 


0.07 


0.25 


2 


2 


> 1 


1 


> 13 


>0.1 


AlO 


9.3 


95 


0.60 


3.2 


0.56 


<60 


>9 


<82 


0.25 


0.07 


0.22 


2 


2 


>0.2 


1 


>2 


>0.02 


BIO 


5.1 


97 


0.65 


10.0 


0.47 


3 


160 


5 


0.25 


0.11 


0.23 


2 


2 


4 


1 


32 


0.3 



" 1 foe =10^^ erg (fifty one erg). 
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Fig. 1. — Orbital separation between black hole and neutron star (solid lines) and specific angular momentum of the gas on the grid (dotted lines) as functions of 
time for Models A2.5, A5, and AlO. The steep drop at the end of the solid lines marks the moment when the neutron star is disrupted. 
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Neutron Star Mass (M^) 



Fig. 2. — Orbital separation as function of neutron star mass for different initial black hole masses and values of parameter a in a simple analytic model (see text). 
Note that the total mass of the system, Mbh + A%s is constant along the lines. Mass transfer leads to orbit widening only for Mbh = 5 and IOMq, whereas GW 
emission decreases the separation. Combining both effects (dashed lines) qualitatively explains the behavior shown for the hydro simulations in Fig. 1. 
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time [ms] 

Fig. 3. — Total neutrino luminosities as functions of time for BH/NS merger Models A2.5, A5, and AlO, and for the NS/NS merger Model A64. 
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annihilation rate, C2.5 




-100 -50 50 100 

d/km 

Fig. 4. — Contours of the logarithm of the azimuthally averaged density distribution (dotted lines) of the accretion torus around the black hole (indicated by the 
white octagonal area at the center) and of the logarithm of the energy deposition rate per cm^ by vv annihilation into e*e^ pairs (sohd lines) for the BH/NS merger 
Model C2.5 at time 9.56ms. The contours are spaced in steps of 0.5 dex. The integral energy deposition rate is 2 x lO'^erg/s. 



